Local topological order inhibits thermal stability in 2D 
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We study the robustness of quantum information stored in the degenerate ground space of a 
local, frustration-free Hamiltonian with commuting terms on a 2D spin lattice. On one hand, a 
macroscopic energy barrier separating the distinct ground states under local transformations would 
protect the information from thermal fluctuations. On the other hand, local topological order would 
shield the ground space from static perturbations. Here we demonstrate that local topological order 
implies a constant energy barrier, thus inhibiting thermal stability. 
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A self-correcting quantum memory [T] is a physical sys- 
tem whose quantum state can be preserved over a long 
period of time without the need for any external interven- 
tion. The archetypical self-correcting classical memory is 
the two-dimensional (2D) Ising ferromagnet. The ground 
state of this system is two-fold degenerate — all-spin up 
and all-spin down — so it can store one bit of informa- 
tion. If the memory is put into contact with a heat bath 
after being initialized in one of these ground states, ther- 
mal fluctuations will lead to the creation of small error 
droplets of inverted spins. The boundary of such droplets 
are domain walls, i.e., one-dimensional excitations whose 
energy is proportional to the droplet perimeter. If the 
temperature is below the critical Curie temperature, the 
Boltzman factor will prevent the creation of macroscopic 
error droplets. Thus, when the system is cooled down 
(either physically or algorithmically) after some macro- 
scopic storage time, it will very likely return to its original 
ground state: the memory is thermally stable. 

This behaviour contrasts with the ID Ising ferromag- 
net whose domain walls are point-like excitations. There- 
fore, they can freely diffuse on the chain at no energy 
cost. As a consequence, arbitrarily large error droplets 
can form, so this ID memory is thermally unstable. 

While the 2D Ising ferromagnet features thermal sta- 
bility, it is vulnerable to static, local perturbations. 
Indeed, an arbitrarily weak magnetic field breaks the 
ground state degeneracy and favours one ground state 
over the other. When this perturbed system is subject to 
thermal fluctuations, the bulk contribution of the mag- 
netic field overwhelms the boundary tension of the do- 
main wall, so once error droplets reach a critical size, 
they rapidly expand to corrupt the memory. This type 
of instability plagues any systems with a local order pa- 
rameter, so they cannot be robust quantum memories. 
Indeed, distinct ground states give different values of this 
order parameter, so a local field coupling to the order pa- 
rameter lifts degeneracy. 

In 2D and higher, there exists quantum systems with 
no local order parameter and whose spectrum is stable 
under weak, local perturbations. These systems have a 
degenerate ground state separated from the other energy 



levels by a constant energy gap, and perturbations only 
alter these features by an exponentially vanishing amount 
as a function of the system size. Kitaev's toric code [2] 
is the best known example of this type. However, ex- 
citations in Kitaev's code are point-like objects — as for 
the ID Ising model — so it does not offer a macroscopic 
energy barrier protection to thermal fluctuations [TJ |3HS] • 
In this work, we study the possibility of combining the 
thermal stability of the 2D Ising model with the spectral 
stability of Kitaev's code to obtain a robust quantum 
memory in 2D. We consider d-level spins located at the 
vertices V of a 2D lattice A = (V,E), with Hamiltonian 

H = -J2 p x, with [P x ,Py] = and \\P X \\ < 1. (1) 

XcV 

We denote the number of spins N = \V\. The term Px 
is supported on the subset X of the spins, i.e., it acts 
trivially on the complement X = V — X of X. The 
Hamiltonian is local in the sense that Px — when- 
ever X has radius larger than some constant w. Since 
we are only interested in the ground state and scaling of 
the energy gap, we can assume without loss of general- 
ity that each Px is a projector. We also assume that H 
is frustration-free, meaning that the ground states min- 
imize the energy of each term of the Hamiltonian sepa- 
rately, i.e., Px\ipo) = IV'o)- Then, the ground space C is 
the image of the code projector P = J\ x Px (henceforth, 
the Px = are not included in such products). 

This family of lattice models, called local commuting 
projector code (LCPC) includes most models of topolog- 
ical order, e.g. [2[71IHj- It has been proved that LCPC 
have a stable spectrum [9HTTj if they obey the following 
local topological order condition. 

Definition 1 (Local topological order). For any topo- 
logically trivial region A, let Pa — Ylx-.xnA^H ^ x ^ e 
the product of projectors that intersect region A. For a 
system with local topological order, the density matrices 
PA = Tr^P and Pa c = Ti^Pa have the same kernel and 
Pa is proportional to a projector. 

Our main result is that any system with local topo- 
logical order has only a constant energy barrier between 
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ground states, similar to Kitaev's code. In fact, our proof 
is inspired by the behaviour of Kitaev's code: we exhibit 
a procedure creating a point-like defect at one boundary 
of the lattice and draging it across to another boundary, 
at a constant energy cost. The difficulty stems from the 
fact that, as noted by Haah and Prcskill [TJ], it is unclear 
whether this process can be realized unitarily in general. 

The rest of the article is organized as follows. We first 
introduce a few definitions and review known facts about 
LCPCs. We then present a local noise model that re- 
quires a constant amount of energy and demonstrate that 
(i) it can be executed in a time that scales linearly with 
system size for any system with local topological order, 
and (ii) it corrupts the information content of the mem- 
ory. We end with a discussion and general conclusions. 

Background — Characterizing the thermal stability of a 
memory requires detailed knowledge of its thermaliza- 
tion process. Since we seek to address a broad class of 
systems, our analysis cannot be model specific. We thus 
retain only two essential features common to all thermal- 
ization processes: (i) the bath interacts locally with the 
system, and (ii) high-energy states are penalized. As we 
now explain, we can combine these features to obtain a 
sufficient condition for thermal stability. 

We will say that a memory with Hamiltonian Eq. ([!]) 
has an energy barrier at most A if there exists a ground 
state if>o and a sequence of T e poly (AT) CPTP maps 
each acting locally on the system and a finite-dimensional 
ancilla A, such that (i) starting from ipo, the sequence 
returns the system to the ground space, (ii) in a state 
that differs from the initial state ipo j an d (hi) the energy 
of any intermediate state is at most A above the ground 
state energy Eq. More formally, these conditions are 

Tr{P-£ T ...£ 2 £ 1 (i> ®pA)] > \ (2) 

Tr{^-£ T ...£2£i^o®PA)]<\ (3) 

max Tr [H ■ £ k . . . £ 2 £ x ® p A )} - E = A (4) 
fcei,2,...,T 

where P is the code projector, i.e. projector onto the 
ground space of H, and the factors | and | are arbi- 
trarily chosen constants. The additional ancillary sys- 
tem, initially in state pa, is used to model some finite 
non-Markovian effects of the bath, so each map £k has 
complete access to it. The energy barrier of a memory 
is taken to be the smallest value of A over all such se- 
quences of maps. If a memory has a macroscopic energy 
barrier A > N a for some constant a > 0, then any short 
sequence of local transformations that returns the sys- 
tem to an altered ground state must visit a high energy 
state, and is therefore thermally stable. Our main result 
is obtained by exhibiting a sequence of maps £k with an 
energy barrier A that is a constant, independent of the 
system size N. In addition, we show that the length T 
of this sequence is proportional to the linear size of the 



lattice when the system has local topological order. 

We call logical operator an operator L that maps the 
ground space to itself, i.e. [L, P] = or equivalently 
[L, Px] = for all X. In particular, we are interested 
in logical operators that act non-trivially on the code 
space, i.e., LP ^ P, as they can alter the encoded infor- 
mation. In a series of paper [T2T - TT5] . it was shown that 
LCPCs always admit at least one non-trivial logical op- 
erator supported only on a ID (constant width) strip of 
the lattice. 

An important subclass of LCPCs are stabilizer codes 
[H], for which Px = ^(I+Sx) with Sx tensor-product of 
Pauli matrices (Jo, 1,2,3 (with o~q being the identity /). Be- 
cause of this particular structure, the non-trivial string- 
like logical operator described above is also a tensor prod- 
uct of Pauli matrices, L — §Z) e k where k labels the £ 
sites along the strip in some natural way, from left to 
right, say. Then, applying the error sequence {cj^} will 
build up to the operator L, and will only visit interme- 
diate states with a constant energy above the ground 
state. Indeed, at an intermediate state < n < £, only 
a segment ®£ of the logical operator L has been ap- 
plied. This segment commutes with all terms Px except 
the ones within distance 0(w) from site k, so only these 
terms contribute to the energy: the excitations are point- 
like objects located in the vicinity on the end of the string 
segment, so the memory is unstable [T31 E3 H3 HE] ■ This 
simple argument fails for more general LCPCs because 
logical operators do not have a tensor product structure. 

Noise model — In this Section we present an error se- 
quence {£k} that achieves a constant energy barrier. To 
motivate our construction, let us first consider a simpli- 
fied version. In a first step, all particles on the strip 
are removed and replaced by particles in random states. 
Mathematically, this action can be described by the max- 
imally depolarizing channel Z\ on all the sites k of the 
strip, which is on average equivalent to applying a Haar- 
random unitary operator on each site k, 

V k [ip] = Tr fe [0] <8 h/D = J U k \i>) (midUk. (5) 

At this point, it is clear (via no-cloning) that the infor- 
mation content of the strip has been wiped out, and that 
it cannot be recovered by any subsequent operation. To 
return the system to its ground space, we can measure 
the projector P, and, with lots of luck, obtain the result 
+1. Of course this outcome is extremely unlikely, and 
for this reason we will refer to this noise model as the 
fortuitous model. Note that prior to the measurement of 
P, we can apply any CPTP map on the strip without al- 
tering our main conclusion that the information has been 
erased. 

Clearly, the fortuitous model is not what we are seek- 
ing. Not only does it rely on an extraordinary amount of 
luck, but the state after the first step has energy propor- 
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FIG. 1: The strip C contains I sites (large circles) whose 
Cartesian coordinates are {(k, 0) 1 < < ^}. Local con- 
straints Pk-i,k act on nearest-neighbours sites fe — 1 and fc 
and on particles in the extended strip C' = \j\ < 1}. 



tional to the strip length. We will now show how to re- 
alize this model sequentially to avoid both of these prob- 
lems. To simplify the presentation, we coarse-grain the 
lattice — i.e., we partition the lattice into balls of radius 
w and view each ball as sites corresponding to a single D- 
level spin with D = d°^ w > — so we can assume without 
lost of generality that: A is a regular £ x £ square lattice; 
the non-zero terms Px in Eq. ([!]) act only on 2 x 2 cells; 
and there exists a non-trivial logical operator supported 
on a single line C of the lattice. Projectors whose support 



intersect L define the strip projector Pc = Yl 
supported on the extended strip C! , see Fig. [I] Similarly, 
projectors whose support intersect sites fe — 1 and fc on 
£ define local constraints Pk-i,k = Ilxn{fc-i fc}^0 ^X- 

The sequential noise model interleaves the depolariz- 
ing and projection steps of the fortuitous model to ob- 
tain individual iterations for every site fc G C, which 
succeed with high probability. Each iteration consists of 
several trials. Trial m of iteration k corresponds to (i) 
applying a trial unitary uj, m ^ on site fc, chosen at ran- 
dom from the Haar measure, and (ii) measuring the local 
constraint Pk-\,k- Trials are repeated until a successful 
trial in which the +1 outcome of Pk-ik is obtained and 
the next iteration begins. Given the state on the strip, a 
unitary is eligible if it leads to a successful trial with non- 
zero probability. The initial iteration k = 1 differs since 
the constraint is not measured. Physically, the whole 
procedure corresponds to creating a random excitation 
at iteration 1, and moving it along the strip across to the 
opposite edge by subsequent iterations. 

The sequential noise model only creates intermediate 
states of constant energy. The reason is the same as for 
stabilizer codes: the excitations are point-like objects. 
Indeed, during iteration k, the state is almost every- 
where indistinguishable from a ground state because it 
obeys all constraints Ps-i.i, except Pk-i, k and Pk,k+i 
since only those potentially do not commute with U^ 71 ^ . 
Furthermore, a failed trial during iteration k does not 
affect the outcome of previous iterations since local con- 
straints commute. Thus, we only need to show that the 
expected number of trials at each iteration is indepen- 
dent of lattice size — so the total duration of the noise 
process grows linearly with the lattice size £ — and that 



the average effect of a complete sequence of £ successful 
iterations has exactly the same effect as the fortuitous 
model. 

Expected number of trials — The sequential model would 
run into a dead-end, an iteration requiring an infinite 
number of trials, if the state of the strip admits no eli- 
gible unitary at the fc th iteration [27]. Such a dead-end 
occurs in the Ising-like toric code introduced in j9] , where 
the plaquette operators B p of the toric code are replaced 
by Ising-like interaction B p B q whose symmetry is broken 
by a single defect plaquette B p * to recover the toric code 
ground state. The sequential model could start prepar- 
ing the B p = — 1 sector and reach a dead-end when it 
encounters the B p * = +1 constraint. However, that code 
does not have local topological order, and its spectrum 
is indeed unstable (B p * is a local order parameter). We 
now show that such dead-ends do not occur with local 
topological order. 

Proposition 2. Local topological order implies that, at 
any iteration k, there exists an eligible unitary Uk- 

Proof. We will prove the contrapositive. Let ip be the 
state during the k th iteration for which no Uk is eligible. 
We have 

Pk-i,kU k \ip) = VU k , (6) 

whose average over the Haar measure, using Eq. ([5|, is 

P k -x,k (Tr fe [ip] ® I k /D) = 0. (7) 

Tracing out the region Rk = {{i,j) : i > k \j\ < 1} C £' 
of the extended strip located at the right of site k, c.f. 
Fig. [TJ Eq. ^ yields 



Tr fe [Pk-i.k] Tr flfc ty] = 0. 



(8) 



Thus, there exists a state |£) in the support of Tr^ [i/>] 
which is in the image of Pi-\^ for i < k but also is in 
the kernel of Tr^ [Pk-i, k ]- This entails violation of local 
topological order on site k — 2 since Tin— i s m the 
kernel of pk-i = Tr^^P but in the image of p 1 fc °f 2 = 

Tr fc=2 [-Pfc-3,fc-2-Pfc-2,fc-l]- □ 

Proposition 3. When the system has local topological 
order, the expected number of trials A k at iteration k is 
finite and independent of the system size. 

Proof. We introduce two maps 

V k -x,k\p] = Pk-i.kpPkM+i (9) 
Qk-i,k\p] = (I - P k -i,k)p(I - Pk-i.k) (10) 

which represent a successful and failed measurement of 
the local constraint Pk-i.k- In a failed trial, the map 
Qk-i.k is always immediately preceded and followed by 
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FIG. 2: Definition of the biasing map E k -\. It is only defined 
on the image of T> k , a projection operation. 



a depolarization of site fc. This sequence can be rewritten 
in an equivalent form, c.f. Fig. [2] 



(11) 



which defines a biasing map £ k -\- This map is not trace- 
preserving since the trace of its unnormalized output 
state is the average probability of a failed trial. 

The sequence of m failed trials followed by a successful 
trial produces the map 

V k - ltk V k {Q k -i,kV k ) m = Vk-i.k {£k-i ® Z?k) (12) 



where we have used Eq. ( 11 ) and V\ — V k . Thus, given 

a state tp, the average probability p k (ip) of a success 

after m failures is p[ m) (tp) = Tr [P k -i,k (^-i ® v k) [ip]] ■ 
Therefore, the expected number of trials 



=i 

[>fc_i, fc ((X fc _i - £ fe -i)~ 2 ® X>fe) [V-]] (14) 



m— 1 

=Tr 



is bounded by the norm of the superoperator inside the 
trace and thus only depends on the microscopic details 
of H, not on system size. Note that the geometric sum 



of Eq. (14) converges since £ k -i cannot have +1 eigen- 
vectors in the groundspace for a local topological ordered 
system, since those would contradict Proposition [2j □ 

Equivalence between models — We now prove the equiva- 
lence of the sequential and fortuitous models. The effect 
of iteration k averages to the map 



Vk-i,k {eiu®v k ) 

rn— 

= Vk-i.k ({i-e k -x)- l ®T> k 



so the total action of the sequential noise model is 

£ = IIfc =2 7 ' fc - i . fc ((^-flfc-ir 1 ® 2 ^) 25 !- 



(15) 
(16) 

(17) 



Terms with non-overlapping support trivially commute. 
We thus move all depolarizing channels to act first, 
which globally depolarizes the strip. To move the bi- 
asing operators past the projectors, it suffices to prove 
that C = [£ k , V k -i >k ] "Dk+i is zero. Because of their non- 
overlapping supports, T>k+i commutes with Vk-\,k and 



C — [£ k T> k+ i,V k -x tk ] — \Pk+lQk,k+l^k+l>"Pk-l,k] — 
since Qk,k+i an d "Pk-i,k commute. Hence, the terms of 



Eq. (17) can be reordered into 



-i,fe 



k=2 



L+1 

IIP 

k=2 



(18) 



k=l 



which has the form of the fortuitous model: complete 
depolarization of the strip, arbitrary transformation on 
the strip, and projection onto the ground space. 

Discussion — While we have focused on systems with 
open boundary conditions, our results extend straightfor- 
wardly to periodic boundaries. The noise process could 
begin at some arbitrary location where it would create 
a pair of defects and, using the techniques we presented, 
wrap one of these defects around the system. A final 
measurement would then attempt to fuse both defects 
back into the ground space. Note that it may be neces- 
sary to wrap one defect several times around the system 
to obtain a successful fusion. For instance, Kitaev's code 
with a twisted boundary requires two complete wraps. 

The noise process we presented corrupts the memory 
after a time that grows proportionally to the system size, 
which can be interpreted as a (macroscopic) upper bound 
to the storage time. However, there are good reasons to 
believe that the actual storage time is in fact independent 
of the system size. At nonzero temperature we expect a 
finite density of defects, so the noise process we described 
could be happening in parallel all over the lattice. As 
pairs of defects meet, they can fuse to the vacuum with 
some probability to create longer error strings. The mem- 
ory time is then related to the percolation of these error 
chains, which should be independent of the system size. 

Our result does not completely close the door to the 
existence of a robust quantum memory in 2D. First, lo- 
cal topological order is a sufficient, but perhaps not nec- 
essary condition for spectral stability. Systems with a 
stable spectrum that do not have local topological or- 
der would escape our conclusions. Second, we have re- 
stricted the form of the Hamiltonian. In realistc phys- 
ical systems, the terms Px need not to commute, the 
ground space can be frustrated, and interaction can de- 
cay algebraically with the distance between sites. Third, 
a macroscopic energy barrier is one mechanism that leads 
to thermal stability, but they may exist other mechanism. 
In particular, a system in contact with a heat bath tends 
to minimize its free energy F = E — TS. Thus, we could 
imagine a system with a large entropy barrier: among all 
possible local noise sequences, only a vanishingly small 
fraction can take one ground state to another, while the 
overwhelming majority lead to dead-ends as described 
above. Such topological spin-glasses Qjjj could offer an 
enhanced quantum memory lifetime. This proposal is 
distinct from existing studies showing that disorder in- 
duces an exponential localization of anyons as 
those only address zero-temperature storage. 
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Conclusion — Our main result hints at a general trade-off 
in 2D between a quantum memory's ability to suppress 
thermal noise and its stability to static perturbations. 
Recent discoveries [2"4Tf2"fj] show that this tradeoff is not 
necessary in 3D. Our result extends prior findings [TBI ITS] 
derived for stabilizer codes to a broader, widely studied 
class of models that includes quantum double [2] , Levin- 
Wen [7], and Turaev-Viro [8] models among others. 
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